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Abstract. The wave functions of a quantum isotropic harmonic oscillator in 
N-space modified by barriers at the coordinate hyperplanes can be expressed 
in terms of certain generalized spherical harmonics. These are associated with 
a product-type weight function on the sphere. Their analysis is carried out 
by means of differential-difference operators. The symmetries of this system 
involve the Weyl group of type B, generated by permutations and changes of 
sign of the coordinates. A new basis for symmetric functions as well as an 
explicit transition matrix to the monomial basis is constructed. This basis 
leads to a basis for invariant spherical harmonics. The determinant of the 
Gram matrix for the basis in the natural inner product over the sphere is 
evaluated, and there is a formula for the evaluation of the basis elements at 
(1,1,. ..,1). When the underlying parameter is specialized to zero, the basis 
consists of ordinary spherical harmonics with cube group symmetry, as used 
for wave functions of electrons in crystals. The harmonic oscillator can also 
be considered as a degenerate interaction-free spin Calogero model. 



1. Introduction 

There are interesting families of potentials invariant under permutations and 
changing of signs of coordinates. The most basic one is a central potential per- 
turbed by a crystal held with cubic symmetry. Another important example is the 
spin Calogero-Moser system of £?-type. In this paper we study a potential which 
can be considered as an iV-dimensional isotropic harmonic oscillator modified by 
barriers at the coordinate hyperplanes, or as a degenerate Calogero-Moser model 
with no interaction. The main object is to study invariant harmonic polynomials, 
which when multiplied by radial Laguerre polynomials provide a complete decom- 
position of the invariant wave functions. This is made possible by use of the author's 
differential-difference operators and the construction of a new basis for symmetric 
functions. 

We consider generalized spherical harmonic polynomials invariant under the ac- 
tion of the hyperoctahedral group Bn acting on R^. Let No = {0, 1,2,...}, the set 
of compositions with N parts is Nq ; for a S Nq let |a| = ^ i=1 on, a\ — Y\.i=i a »- 
and (t) a = l\ l=1 {t) a ., (the Pochhammer symbol is (t) n = n"=i (* + * - !))• De- 
note the cardinality of a finite set E by #E. The set V of partitions consists 
of finite sequences A = (Ai,A2, •••) € Ujv=i satisfying Ai > A2 > . . . > 
and two partitions with the same nonzero components are identified. Then 
V n = {A e V : |A| = n} for n = 1, 2, 3, ... . For a partition the length / (A) — 
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# {i : Xi > 1} is the number of nonzero components. For a given N we will need 
the partitions with I (A) < N, denoted by V {N) = VnNtf , and V { n N) = 7>„nN^. 
One partial ordering for partitions is given by containment of Ferrers diagrams: 

thus fj, C A means /i, < A-j for each i. The notation = ( 0, . . . , 1, . . . ] for the i th 



standard basis vector in Nq (for 1 < i < N) is convenient for describing contiguous 
partitions; for example A + E\ = (Ai + 1, A2, • ■ •)• 

For x G R N and a G let x a = TT^Ii , a monomial of degree \a\. Let = 
span {x a : a G Nq , \a\ = n}, the space of homogeneous polynomials of degree n > 
in N variables. The Laplacian is A = J2iLi an d the Euclidean norm is 
■A 1/2 

xf . Fix the parameter n > 0. 



We will be concerned with the operator A + 2K^ i=1 t^rgf-- It is associated 
with the Calogero-Sutherland model with the potential function V (x) = ui 2 \ \x\\ 2 + 
N 1 

k (k — 1) ^2 ~ ( w ith u) > 0). This is an iV-dimensional isotropic harmonic oscil- 

»=i %i 

lator modified by barriers at the coordinate hyperplanes; or it could be considered 
as a degenerate form of the type-i? spin model for N particles with no interaction. 
The spin model with interactions and reflecting barriers was studied by Yamamoto 
and Tsuchiya O]. With the base state 



N 

-1 

we obtain the conjugate of the Hamiltonian, for any smooth function / on 



V(x)=exp (-|||x|| 2 )n 



N 1 \ 

i=i * i 



^(xY x -A + c 2 ||x|| 2 + K ( K -l)> ) (*)/(*)) 



= ^-^E-^ + ^^ + ^ + ^E^ /(*)• 

\ i=l 1 1 i=l l / 

The version of the differential-difference operators introduced by the author 
(||, or see || Ch. 4]) are the fundamental tools for analyzing this eigenfunction 
problem. The action of the group Bn on M. N induces an action on functions, 
denoted by wf (x) — f (xw) for w G £>at. For 1 < i < N let <7j denote the reflection 

(xi, . . . ,x N ) a l = (xi, . . . , —Xi, ■ ■ ■ ,x N ) , 
and define the first-order operator T>i by 

for sufficiently smooth functions / on M. N . The Laplacian operator is 

JV 

Ak=E^' 

i=l 

then 

A K f (x) = Af{x) + n}^\-—f{x) - 2 

i— 1 ^ 1 1 % 
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If / is even in each xi (that is, -invariant) then A K f (x) = Af (x) 
2nJ2iLi ~t~ ~sk~f ( x )- ^ we modify the above Hamiltonian to 



H= -A + oj 2 \\x\\ 2 + K Y, 



N 

K - 



which is identical to the previous one when acting on -invariant functions, we 
obtain 



>■■<■}{> U i = ( - -A,, + to^Xi-^+wN (2k + l))f(x). 



The eigenfunctions of this operator can be expressed as products of radial Laguerre 
polynomials and homogeneous harmonic polynomials. 

Definition 1. For n > the space of harmonic (sometimes called h-harmonic) 

homogeneous polynomials of degree n is H n = |/ G : A K f = o\ . Let = 

{/ £ Hbn : o"j/ = /, 1 < i < iV}; i/izs is i/ie space ofTL^ -invariant harmonic polyno- 
mials, even in each Xi,l <i < iV. 77ie subspace o/H^ consisting of polynomials 
invariant under permutation of coordinates is denoted by . 

It is not hard to show that if / € H n then i^ 1+n / 2+Nk x ) ^ w ||x|| 2 ^ / (x) is an 

eigenfunction oi_ip~ 1 Hip with eigenvalue a; (2n + 4s + N (2k + 1)) (this calculation 
uses equation 3.1). So combining bases for each H n with Laguerre polynomials 
provides a basis for the polynomial eigenfunctions (positive energy) of vp~ 1 'Hip- 
The use of functions of this type in the general Calogero-Moser models of types A 
and B is discussed in van Diejen |Q. 

The classical motion problem corresponding to TL is easily solved: in the one- 
dimensional case the particle with mass 1 at s € K satisfies (4^ s (t) — — J^V (s) = 
— 2w 2 s + 2k(k — 1) s~ 3 with the solution (for s > 0, k > 1) 

s(t) = (q + a sin (2 3/2 uj (* - *o))) ^ 



q= (a 2 + k (k - 1) /uo 2 ) 



1/2 



here to is an arbitrary phase shift, a > is arbitrary, and the energy is | + 
^(s) = 2 9 w 2 . 

Primarily we will study the invariant functions associated with A K ; the invari- 
ance is with respect to the Weyl group of type B, thus the functions are invariant 
under sign changes {oi\ and permutations of coordinates. Such functions are ex- 
pressed as symmetric functions of the variable x 2 — (x 2 , . . . , Xjv). Note that a 
special case of this study is the problem of spherical harmonics on R 3 which sat- 
isfy i?3-invariance, appearing in the wave functions of electrons in crystals. We 
will construct an explicit basis for these polynomials, but unfortunately it is not 
orthogonal. There seems to be a good reason why an orthogonal basis has not yet 
been found: let i?y = x{Dj — XjT>i (then RijA K — A K Rij, and Rij is an angu- 
lar momentum operator), and let S n = X)i<j<j<Ar ^Hp) then Si is the Casimir (or 
Laplacc-Beltrami) operator, for each n the operator S n is self-adjoint in the natural 
inner product on the sphere (defined later) but £2,53 do not commute, already for 
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./V = 3. So the usual machinery for constructing good orthogonal bases (like Jack 
polynomials) does not work here. The conjugate 

Xi Xq \ 

^R.^ = x x K I a . _ ai \ 

commutes with the Hamiltonian Ti (see Taniguchi || for a more general treatment 
of commuting operators in the context of r~ 2 type potentials). By straightforward 
calculation 

Si = R ij 

l<i<j<N 



Wl 2 



Thus Si has the eigenvalue -2n (N - 2 + 2n + 2Nk) on E% n . 

2. A BASIS FOR SYMMETRIC FUNCTIONS 

In this section the number N of variables is not specified, with the understanding 
that it is not less than the length of any partition that appears. We use the notation 
of Macdonald Q for symmetric polynomials in the variable x G M. N . Let Sn denote 
the symmetric group on N objects, considered as the group of N x N permutation 
matrices, acting on the left on Nq . For any a € Nq let a + denote the unique 
partition wot G "P^, for some w € Sn (the sorting of the components of a in 
nonincreasing order; w need not be unique). 

Definition 2. For A G V^ N ' the monomial symmetric function is 

m x = J2{x a : a G f<\a+ = A} , 

summing over all distinct permutations of A. The elementary symmetric function 
of degree 1 is 

N 

ei = J2xi. 

i=l 

It turns out that the basis elements for invariant harmonics are labeled by par- 
titions with the property Ai = A2; that is, dimH^ = # |a G T , h N ' > : Ai = A 2 |. 

Further the formula for the projection onto harmonics uses powers of ||x|| 2 = 
&i {x\, . . . ,x N ^j . This leads to the following definition of a basis for symmetric 
functions well suited for the present study. 

Definition 3. For A G the modified monomial symmetric function is 

"^a = ei 1_A2 TO (A2 , A2iA3i ... ) . 

We will show that {m A '■ A G is a basis for the symmetric polynomials 

on WL N and the transition matrix to the m-basis has entries in Z, is unimodular 
and triangular in a certain ordering. One direction is an easy consequence of the 
multinomial theorem. We use the notation ( J J for the multinomial coefficient, 
where a G 7* N , J2iLi a i = 3i an d („) = if each ai > 0, else (•') = 0. 
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Definition 4. For A, v £ £/ie coefficient (*\ is defined by the expansion 



|A|=M+J 



A^ 



/or j = 0,1, 2,. 



We establish some basic properties as well as an explicit formula for this coeffi- 
cient, and also explain the relation with the generalized binomial coefficient. 

Proposition 1. For A, v g the following hold: 

1. /ei j = |A| - > 0, then 

note (\£ ) = if any component Xi — o~.- L < 0, the sum is over all distinct 
permutations of v, (including any possible zero components to make v an 
element of Nq ); 

2. ( A ) ^ Q if and only if v C A; 

3- 6 = 1- 
Proof. Expanding ejm„ we obtain 

4 mu = { Q ^ : a € N^, |a| = j, ^ 

Since this is a symmetric polynomial, the coefficient of m\ equals the coefficient of 
x x , namely the sum of (^) with a + a = A (component- wise addition) . For any a 
with a + = v \i Oi > Xi for some i then the corresponding term is zero, by definition 
of the multinomial coefficient. For part (2), v C A (and |A| = \v\ + j) implies 
> 0; conversely if ( A f_ cr ) ^ for some a with a + = v then v C A (if o~j < A; 
for each i then cr + C A; indeed, let w be a permutation so that vi — cr w u> for each 
i, so that Vj <Vi = u w u\ < X w u-\ for 1 < i < j then Uj < Xj because it is less than 
or equal to at least j components of the partition A). Part (3) is trivial. □ 

The term "generalized binomial coefficient" has appeared in several contexts 
(such as Jack polynomials, see Lassalle @] ; part 2 of Proposition [l] is a special case 
of his results); here we use the version defined by 

m x (x + sl N ) ^ m v( x ) 



(2- 1 ) v „ NS ' = > s w - M [ v ;,' s£R. 

m\(l N ) \v)m v {l N y 

This is essentially the Taylor series for m\ evaluated at (xi + s, . . . ,%n + s). The 
coefficient [ J = \ V This formula is established by means of the 



y isJ vl(\\\-\u\)l\v 

Let (f,g) = 
and set s = to show 



inner product (f,g) = f (d) g (x) \ x =ol a Pply iw'Y to both sides of equation 2.1 
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take the inner product with ra a for any a G with |er| = |A| — j; note that 

(m a ,m u ) = 5 av a\m a (l N ) and ^7) ™a,«V^ = ^m^e^m^. The coef- 

ficient ( A ) is thus shown to be independent of iV, provided the number of (nonzero) 
parts of A, v does not exceed N. There are useful formulae involving contiguous 
partitions (if v D A and \v\ = |A| + 1 then v = A + Si for i = 1 or A;_i > A;). 

Proposition 2. Suppose X,a E V and \a\ > |A| + 1, also j = 1 or Aj_i > Aj iften 
(2.2) ^ A + £j j = l + #{ 4 :A ! = A J + l}, 

m (:)=E{( A ; ei )( A r>--^>*< 

Proof. Suppose A s _i > A s = Aj_i = Xj + 1 (or ,s = l,or j = 1), then there are 
exactly j — s + 1 distinct permutations wA for which A + Sj — wX has no negative 
components (namely the transpositions w = (i, j) for some i with s < i < j — 1, or 
w = 1). Each such term contributes 1 to the sum defining ( ~t e ')- For the second 
part, let n = |o~| — |A| and consider the coefficient of m a in 

AY™ _„»-iv^ + £ 

|i/|=|A|+n 

x -> x / V \ /A + £ 



e 1 m x = 2^ \X/ mu=ei 2s\ X 



|i/|=|A|+n » 

summing over z = 1 and i > 2 with Aj_i > A,. □ 

We are ready to consider the transition matrices between the m\ and rh\ bases 
(that the second set is a basis will be immediately proven). 

Definition 5. The transition matrices A((i, A) , B (ji, A) (for X, fi £ T 3 , |A| = 

are defined by 

fh\ = ^A(fj,,X) m^, 
m\ = B (fi, A) m M . 
By Definition [I| it follows that 

^ ' ^ ^(^2, A2, A3, 
The condition A (/i, A) ^ leads to the following ordering of partitions: 

Definition 6. For A,/x G P i/ie relation fi -< X means that \X\ = |/Lt| and Xi < /i; 
/or all i >2. 

It is clear that ^ is a partial ordering, /1 ^ A implies /ii < Ai , Z (A) < I (/i) and 
that the maximum element of is (rt). In fact /i ^ A implies that A domi- 

nates /1, but ^ is not identical to the dominance ordering. Further the partitions 
{A : Ai = A2} are minimal elements; indeed suppose /1 ^ A and Ai = A2, then 
Mi < Ai = A2 < M2 which implies fx± = {12 = Ai; thus J2i>2 (/■*« ~ ^i) = and 
X = fi (by definition, /ii — Ai > for each j > 2). There are exceptional minimal 
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elements of the form (to + 1, m, . . . , m), with n = Nm + 1, but these are minimal 
only for since (to, to, . . . , to, 1) ^< (to + 1, to, . . . , to, 0). We see that the ma- 

trix A (/i, A) is triangular (A (/i, A) ^ implies // r< A), unipotent (A (A, A) = 1) and 
its entries are nonnegative integers. 

Proposition 3. For n > the set is a basis for the symmetric 

polynomials of degree n in N variables. The transition matrix B (fi, A) is triangular 
for the ordering ^ and is unipotent with integer entries. 

Regarding the dependence on N: begin by assuming that n < N so that I (A) < 
N for each A <E V n and the defining equations for A (/x, A) , B (/x, A) are unambiguous, 
then to restrict to a smaller number of variables, say M < N, substitute Xm+i = 
xm+2 = ■ ■ ■ = Xn = with the effect of removing all terms with I (A) > M or 
I (fx) > M. Thus the transition matrices for the case of M variables with n > M are 
principal submatrices of the general transition matrices (deleting rows and columns 

labeled by A with I (A) > M). If the index set Vn is ordered first in decreasing 
order of I (A) and then with respect to -< (possible, since fx ^ A implies I (fi) > I (A)) 
then the triangularity makes it obvious that the submatrix of B is the inverse of 
the submatrix of A; see the example of B for V\ at the end of this section. 

It remains to establish an explicit formula for B (fi, A). The formula itself was 
postulated based on computer algebra experimentation, and explicit known formu- 
lae for N = 3. The proof will be by induction on n and requires showing that the 
claimed formula satisfies the following contiguity relations. 

Proposition 4. For A, fx S V n , and a S V^X with o~\ = 02, 
(2.4) Yl rW + Ei,f)=B0*.A) l 

iOA,M=|A|+l 



(2.5) £ r)B(a,u)=0. 

iOA,|i/|=|A|+1 * ' 

Proof. Multiply both sides of the equation m\ = J2n-<\ B (fx, A) to m by e\ to obtain 



iOA,|i/|=|A|+1 aiXA 



Y {\ ) mv = ^2 B (t 1 ' A ) m »+zi 



iOA,|i/| = |A| + 1 x 



Since the set |to t : t G is a basis the desired equations are consequences of 

matching coefficients in the two right-hand expansions. □ 



Lemma 1. Suppose B' (fx, A) is a matrix satisfying the equations 2.4,2.1 with B 
replaced by B' , has the same triangularity property as B,and B' (a, a) = 1 for each 
creP with o\ = (72, then B' (fx, A) = B (fx, A) for all /i,AeP (with |A| = \/x\). 

Proof. This is a double induction on |A| and Ai — A2. For any a € V with o~i = o~% the 
hypothesis shows that B' (fx, a) = for fx ^ a (since a is minimal) and B' (a, a) = 1, 
thus .B' (/i, a) = B (/i, cr) for all (with = |oj). Suppose that i?' A) = i? (/z, A) 
for all fx, \ £ V with |A| = = n or with |A| = \fx\ = n + 1 and Ai — A2 < j for 
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some j > 0. Fix v G V with \ v\ — n+l and v\ — v% = j + 1 and let A = z/ — £i. Then 
r D A and |r| = n + 1 implies r = v or r = A + Ei with i > 2 and Aj_i > A; = i/,-. 
Note that ( A ^ £1 ) = 1 by equation |2.2| . For any /i £ "Pn+ij by hypothesis 

'X + ei 
X 



E 

i>2,A i _i>A i 



B' (ji, X + Ei), 



replacing B' (/x — e%, A) by if (i\ = (i2- By the inductive hypothesis each term 
B' (<t,t) appearing on the right hand side satisfies B' (a, r) = B(<j,t). By the 
Proposition, B' (fi, v) — B (/i, v). This completes the induction. □ 

In the following formula the coefficient (Z) is used with partitions whose first 
parts are deleted, also (/12 — I, ^3, ... ) is not a partition if (12 = M3 and so the a + 
notation is used; the binomial coefficient ( _1 ) = 0. 



Theorem 1. For (i, X e V n 



B(/i,A) = (-l) 



Ai — /ii 




Ai - fi 2 - 1 



(M2 - 1,M3) 
(A2, A3, . 



Proof. The proof is broken down into cases depending on the values of Ai — (i\ 
and (i\ — H2 ■ We let B' (//, A) denote the right-hand-side of the stated formula and 
proceed as in the Lemma to show B' satisfies equations 2.4 and |2.5| . It is clear 
that B' (/i, A) 7^ implies < fii for all i > 2, that is, jit ^ A. Also B' (/i, (i) = 1 
for each ^ g TV Fix (i, X e V n and let // = (/x 2 , /i3, ■ ■ • ) , m" = (M2 — 1, M3> • • • ) , 
<r' = ((T2, 03) ■ ■ ■ ) , <j" = {<72 — 1, 03, . . . ) + , A' = (A2, A3, . . . ) (these are partitions, 
but with components labeled by i > 2; for example A' +e 2 = (A2 + 1, A3, ...))• The 
partitions v satisfying i/DA an d \v \ = |A| + 1 are of the form A + £, with i = 1 or 
Ai-i > A^ We rewrite equation 2.4 (not yet proven!) as 

'X + Si 



i>2,Ai_i>Ai 



A 



B' ((i + £1, A + ei) = B' (p, A) -B'(n + £1, A + £1) . 



The left-hand-side equals 



Al - (12 

v /il + 1 - /i 2 

The right-hand-side equals 



i>2,A i _ 1 >A i 



Ai 
Mi 



M2 
1 - 



- 1 
M2 



A* 

A' + Ei 



A + £, 

A 



-1 



= (-1 



Ai 
Mi 

Ai - M2 - 1 
Mi - M2 

Ai— AH— 1 



Ma\ _ Ai + 1 - Ma 
M2/ \Mi + 1 - M2 

Al - (12 

Mi + 1 - M2 

Al - (12 
Ml + 1 - ^ 2 



Ai - M2 - 1 
Mi + 1 - /i2 



Similarly rewrite equation 2J5 as 

/A + Ei 



i>2,Ai_i>Ai 



A 



S' fa-, A + Ei 



-B' (tr,A + ei). 
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with the left-hand-side (of course, <j\ = o~ 2 ) 
(-l) Al ~ CTl 

'Ai — (TaX / a' \ /Ai — <j 2 — 1\ / a" \ \ /A + e, 



x E 

i>2,Ai_i>Ai 



K a\ — o<x) \A' + £,/ V °i - CT 2 /\A' + et/J\ A 
and the right-hand-side 



■1) 



Ai -cti 



Ai - a 2 \ /a'\ A - <r 2 - 1\ la' 
ui-a 2 )\\') \ ax-G 2 A A' 



The goal is to show that the two sides are equal. We reduce the two cases to one 
by replacing ai,a',a" by /ii + 1, //, //'respectively. Let i? denote the difference 
between the two sides (left - right), and let 

{ \l-(i<2 \ /n' 



A*i + 1 - M2/ \A 

S= A Ai-/x 2 \ In' \ /A + sA _ A-/x 2 -1Vm" 

1 VMi + l-W .> 2 ^ i>A AA' + ei/\ A / U1 + I-M2 



£ 2 = 



Ai — M2 — 1\ / M" \ / ^ + £ 



' ii + 1 -^A> 2 ,Ar 1 >A 1 \ A, w\ A 



then £7 = (— l) Al Ml 1 + J3i + E 2 ). We can assume /«' D A' (// ^ A), because if 
any term in E is nonzero then (^,) ^ (if ( A /^ £ .) 7^ for some i then \j! D fi" D 
A' + Si D A', and similarly if ( A /^ e .) 7^ for some i then // D A'). 

1. Case: |//| = |A'|: Here E\ = = £?2, because any coefficient (^) with 
|cr| < |t| is zero. Further = n — |//| = Ai and thus i?o = (from the 
binomial coefficient). 

2. Case: = |A'| + 1 and /ii > /Lt 2 : Here £2 = and Ai = fii + l,thus 
the coefficient of (^") is zero (specifically, ( Xl ^ 2 ~ 1 ) - ( 2 ) = 1 - 
1). The hypothesis implies fi' = A' + £j for some j > 2, that is, /i = 
(Ai - 1, A 2 , . . . , Aj ■ + 1, . . . ). If j = 2 then Ai > A 2 + 2 otherwise j > 2 
and Ai - 1 > Aj_i > Aj. Thus E + Ex = ( X+ X £j ) - C^ 3 )] only the term 
with i = j in the sum is nonzero (we showed Aj_i > Aj). By equation 2.1 

< A+ / J ) = ( X 'P), since 1 <£ {i : A* = A, + 1}. 

3. Case: |//| = |A'| + 1 and /ii = /i 2 — 1 (that is, \i — a — £1 where ci = 02): 
Here £?2 = and the binomial coefficients in Eq,Ei are all equal to 1. Let 
\j! = A' + Ej. If fi" = A' then either j = 2, A = (Ai, Ai — 1, A3, . . . ) , fx = 
(Ai - 1, Ai, A 3 , . . . ) with Ai > A 3 or j > 2, A = (Ai, A x , ... , Ai, Ai - 1, Aj +i , . . . ) , 
fi = (Ai — 1, Ai, . . . , Ai, A J+ i, . . . ) (that is, A^ = Ai for 1 < i < j — 1 and 
Xj = Ai — 1), which implies 



Eo + Ei — 



fi' \/\ + £j\ I n"\ j (if 



A' + £j/\ A / \A7 \A' 



by equation ^2[ If (i" ^ A' then Aj ■ + 1 < A 2 (so j > 2) and E + Ei = 

( A+ A^)-( A »0. 
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4. Case: > |A'| + 2 and /xi > (i 2 : Thus Ai > p,\ + 2 > fi 2 + 2 > A 2 + 2 and 
equation 2.2 implies ( "^ g * ) = ( ) for each i > 2, Ai-i > Ai (which includes 
i = 2). Also equation |2.3| applied to the truncated partition t' = (r 2 , T3, . . . ) 
with |r'| > I A' I shows that 

Substituting r' = // and r' = /i" in this equation shows that E0+E1+E2 = 0. 

5. Case: > |A'| + 2 and ^1 = ^2 — 1 (that is, /x = a — e\ where <j\ — <r 2 ): 
Similarly X 1 > (fi 2 - 1) + 2 > A 2 + 1, and ( x \ £ ') = ( X 'y €l ) for each i > 2 
with Aj_i > A 4 . Either A x > A 2 + 2 which implies ( A ^ £2 ) = ( A 'J £2 ) = 1 
or Ai — A 2 + 1 which implies /i 2 = A 2 and ( A /^_) = = \^{\ e ) because 
(A' + e 2 ) 2 > /i 2 > (/i") 2 - ln both cases equation |7| shows i?o + E\ + i? 2 = 
as previously (in the case \i = A 2 + 1 one has ( X \ £2 ) = 2, ( A J £2 ) = 1). 

Thus the matrix B' (A, /x) satisfies the hypotheses of the Lemma, and the proof 
is completed. □ 

Case (5): |//| = |A'| + = a — E\ with o\ — er 2 , is the only case in the 
proof which depends on the \i" part of the formula (of course the first steps of an 
induction proof are crucial!). By way of an example, here is the matrix B for Vq 



1 


-2 





-2 


8 





2 


-18 


18 





1 








-10 








42 


-72 








1 


-3 


3 





3 


-9 


9 











1 


-3 





-2 


13 


-18 














1 








-9 


24 

















1 


-2 


2 


-2 




















1 


-4 


9 























1 


-6 


























1 



with the rows and columns labeled by [2211,3111,222,321,411,33,42,51,6] (sup- 
pressing the commas in the partitions). 



3. Invariant harmonic polynomials 

In the sequel the symmetric polynomials have the argument x 2 = (x\, . . . ,x^) 
(and y 2 is similarly defined). Thus e\ (a; 2 ) = ||x|| 2 , and the former will often be 
used in equations involving symmetric functions. We recall the basic facts about 
the Poisson kernel, that is, the reproducing kernel, for harmonic polynomials H„, 
(see Dunkl and Xu Ch. 5]). The intertwining operator V is a linear isomorphism 
on polynomials such that VP^ = for each n, T>iV = V-^r- for 1 < i < N, 
and VI = 1. Here the reflection group is a direct product and so the intertwining 
map is the iV-fold tensor product of the one-dimensional transform (see || Theorem 
5.1]), defined by (for n > 0) 

Vx\ n — x 2n y x ^ n + 1 — Wn+l Jln+\ 



Tt + 1 
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Then let K n (x,y) = ((x,y) n ), for x,y G R N ,n > (V^ acts on x). The 

key properties of K n needed here are K n (xw,yw) = K n (x,y) for w G -B/v and 
£>f (a;, y) = y t K n _ 1 (x, y) for 1 < i < N. 

Definition 7. Let du> denote the normalized rotation-invariant surface measure on 
the unit sphere S={i£ M. N : \ \x\ \ = l}, and for polynomials f,g the inner product 
is 

N 



(/,s)s = c « / f (x) g (x)Y[\x l \ 2K da; (x) 



where c K = — — — -^jf- so that (1, l) g = 1. 



r(#)i>, _,, 



The reproducing kernel for H„, n > 1 is 



'iV 

0, 2/) = ( TT + Nk + n ~ 1 



X 

j<n/2 



Nk] 

£ V2 ' ' J *- 2J INI 2 ' IM| 2i ^- 2j (*,y) 



and satisfies (P n (-,y),f) s — f (y) for each / G H„. The general theory g] for 
our differential-difference operators shows that / G H n ,<7 G H m ,m 7^ n implies 
(/>5)s = 0- Let denote the kernel i^2n (x,y) symmetrized with respect to the 
group 1*2 (f° r nxe d x G M. N sum over the 2 N points (±Xi, . . . , ±zjv) and divide 
by 2 W ); it is the reproducing kernel for H^. Let K® n (x, y) , P^ (x, y) denote the 
symmetrizations of Km {x, y) , Pm (x, y) respectively, with respect to the group B^. 
Thus Kg (x, y) = £t J2 wgSn K$ n (xw, y). 



Proposition 5. For n > 1, x,y G 



piV 



K 2 ° n (», „) = 2- 2 « £ „, r , n * 2 V° 



aGN",|a|=n " v '" ' 2 '" 

K ^ V) - 2 A! ' 

Proof. By the multinomial theorem, =^yy E { { 2 p) x0 V = 2n|. 

Symmetrizing over removes monomials with odd exponents, and thus -f^n (- t j 2/) 
is the result of applying V to J2 {^J^y x2a y 2a '■ a e ^01 M =n\- Further (2a)! — 
2 2 H a | (|) a , and ^i 2a = ((|) q / (« + |)J x 2q . For a fixed A G the sum 

M E {x 2wa y 2a : u> G Sjy, a + = A} is a multiple of m\ (x 2 ) mj, (y 2 ), and evaluated 
at x = l N — y the sum equals m\ (l N ) (that is, |{a£ Nq : a + = A}). The terms 
a! and (k + 4) are invariant under Sjy. □ 
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Corollary 1. For n > 1, x, y e 



J>£ (*, y) = (y + JV« + 2n - l) £ U+ ? , ^ i (-T ex (y 2 ) 

E 



euu; ; ei V y 2V 

3=0 

1 m x (x 2 ) m A (y 2 ) 

(JV) 



\ ;_A!(.+ I) A m A (l") 



Let Pf„ denote the space of £?Ar-invariant elements of then dimP- 8 — 

(N) 



- 2n 

#Vn" ' ■ The space H^ n of -Bjv-invariant harmonic polynomials is the kernel of A K ; 



since A K commutes with the action of B N and maps P 2 ^p onto P 2 ^l 2 we see that 
dimH^ = #vi N ^ — #V^!\. The map A A + e\ is a one-to-one correspondence 
of V^li onto a subset of whose complement is 

pW = {AepW:A 1= A 2 }. 

Thus dimHf„ = #vi N \ and we will construct a basis for whose elements are 
labeled in a natural way by vi N ^ . There is a generating function for the dimensions: 
for A e Vn N ^ the conjugate A T (the partition corresponding to the transpose of the 
Ferrers diagram) of A is a partition with 2 < (A T ) . < N for all i (no parts equal to 
1 or exceeding N); thus 

oo N 

E(#^ v) )* B =ii( i -*r 1 - 

n=0 j=2 

This expression yields an estimate for ffP^. Indeed, let M = 1cm (2, 3, . . . , N), 
11^2 (! - 1 = P (?) (! - g M )~ (jV_1) for some polynomial p (q); thus #p( N) = 
O (^(jj) N 2 ) as n — > oo. In (x, y) expand each m A (j/ 2 ) in the m^-basis 

K 2n {x,y) 2 ^ ^ ^ + ^ ^ (iJV) m„ (y ) , 
this leads to the following. 

Definition 8. For ii G vL N ^ the Bn -invariant polynomial is given by 



(JV) 



Recall that A e P„ and £> (A, /x) 7^ implies I (A) < I (/u) and AeP, 

Theorem 2. For /1 G Vn N \ if Hi > M2 A re /i^ = 2/i^_ £l and i/ = /li 2 
A K /i M = 0. 

Proof. From the basic properties of if 2n it follows that A« x ^ 2 „ (x, y) = | |y 1 1 2 K 2n -2 ( 
Symmetrize this equation with respect to Sjy (and note that A K commutes with 
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the group action) to obtain 



Men 

,.2\ isB 



= ei (y 2 )Ki n _ 2 ( x , y ) 
= 2 l - n J2 K(x)m a+Sl (y 2 ). 

By definition e\ (y 2 ) fh a (j/ 2 ) = m CT+£l (y 2 ) . Considering the equations as expan- 
sions in |m M (y 2 ) : /i € T 3 ^'' 1 j shows A K h a+£l (x) = 2h a (x) for a 6 "P^i and 
A K /i M (x) = if fxi = /ia- □ 

Corollary 2. TTie set : fj, E vi N) } is a basis for . 

Proof. For X, fx E V { n N) the coefficient of m A in is 5 AjU (2"A! (k + |) A m A (l^))" 1 
and thus : jj, E "Pn^H is linearly independent. Also dimHf n = jfiP^K □ 

Besides the inner product on polynomials defined by integration over the sphere 
S there is also the important pairing defined in an algebraic manner using the 
operators T>i , namely 

{f,9) h = f(V 1 ,...,V N )g(x)\ x=0 . 

Since Vfx 2n = 2n(2n- l + 2n)x 2n ~ 2 and V 2 xf L = for j ^ i, we have that 
(x 2n ,xf l ) h = 2 2n n\ (k+ |) n and (x 2a , x^) h = 5 ct/3 2 2 l a l a! (/s+|) a for a,/? £ Nff, 
so that monomials are mutually orthogonal. It follows that (rn\ (x 2 ) , m M (^ 2 )) ;i = 

<5 Ai u2 2n A! (« + |) A m A (1^) for A,/j£ P^. It was shown in ( |§ Theorem 3.8], see 
also §, Theorem 5.2.4]) that for f,g E H n 

This allows the direct calculation of (h\,h^) s for A,/i G 
Proposition 6. for A, /i E Vn 



Ih h \ ST g_(A 1 ojB_(/x 1 oO 



(/ iA ,/ lM ) s = 2- 2 " ^^+JVk) ) 

Computations with small degrees (n < 8) show that there are no product (linear 
factors in n) formulae for the inner products and {h\} is not an orthogonal basis. 
It is however possible to find a nice formula for the determinant of the Gram matrix 
[(h\, hf^fj A ^ e p{N). Before presenting this result we consider the bi-orthogonal set 

for ^hp : fj, € vL N ^ constructed by extracting a certain multiple of the coefficient 
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of m\ (y 2 ) for A S in the Poisson kernel P 2 ^ (x, y). Indeed P® n is a multiple 

of 



E E (-DV^ei^ 



J=0 „ 67 ,< 



= E E E 



•>>i ( f + NK ) 2n -\-j (x 2 ) m a (y 2 ) 
x 2\3 I A\ (f + M^n-i-j m CT (x 2 ) m A (y 2 ) 



a/ («+§)„ m CT (l w ) 



Then let 

_ A ei (x 2 ) J _ /A\ 1 m g (x 2 ) 

ffA(X) (_A_jVK-2n + 2). £f W*!(«+|) 

j=0 ■> \ 2 /j ff£p (f() ffCA » ' V 2/cr 

so that the coefficient of m>, (y 2 ) in -P^j (^i u) is (y + ^ K )2n 9 X Observe that 

9x ^ = \\ ( K + 1) mx (pV) + 1 1^1 1 2 5a l^) where g A (x) is a polynomial of degree 
2n- 2. 

Proposition 7. For A, /i S Pri^' inner product 

(5X '^ )ft = ^A! ( K +I) A m A (l^)- 
Pnoo/. Indeed ( flA , h^) h = ^—^—^ ( mx ( x 2 ) , h„ (x)) h + g' x (V) A K h„ (x) = 
^ A!(K+^)"m A (i w ) ' by thc P ro P crtics 01 thc Pacing (;-) h - n 

Essentially g\ is (a scalar multiple of) the projection of m\ (x 2 ) onto H2„. The 
method of orthogonal projections to construct bases of harmonic polynomials was 
studied in more detail by Y. Xu JTofl . 

We turn to the evaluation of the determinant of the Gram matrix of {h\ : A G 

Vn^} using the (•, inner product; the value for (•, -} s is just a product of a 

power (#Vn ) of the proportionality factor with the previous one. It is perhaps 
surprising that the calculation can be carried out by finding the determinant for 

the entire set , which can easily be done. By the orthogonality of 

{m^} we see that the Gram matrix G (A, fx) = (h\, h fl ) h for A, £ Vn satisfies 

G = BCB T , 
det G = (det B) 2 det C = det C, 

where the diagonal matrix C (A, fi) = Sx^ (XI m x(^ N )) \ just as in 

Proposition | (h\,h tJ .) h = T^JY^Tm for an y S vk N) ■ We al- 

ready showed B is triangular with 1 's on the main diagonal so that det B = 1 . 
Let 



D n = det C = 

xev. 



n>HW"))' 
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and let D n = det G, where G is the principal submatrix of G for the labels A, [i S 
Vn . We will show that 

D n = ^- J] ((A 1 -A 2 + l)(AT/«+y + 2(n-l)-(A 1 -A 2 ) 

A simplified form of this will be given later. The idea of the proof is to use the 
orthogonal decomposition P„ = ©j<„/ 2 1 1 x 1 1 2j ^-n-2j to produce a transformation 
of the Gram matrix into block form. The underlying relation is the product formula 
for A K : for to, j > and any / € Pm 

A K |M| 2J f( x )=Aj(m + j-l + ^ + NK^ \\x\\ 2j - 2 f(x) + \\x\\ 2 "A K f(x). 

Specializing to harmonic polynomials / £ H m and iterating this formula shows that 

(3.1) A* | \x\ \ 2j f (x) = 2 2 * (-i) s (-m - j + 1 - y - | |x| \ 2j ~ 2s f (x) 

for s > 0; note A 6 K | \x\ \° f (x) = for s > j (see ||, Theorem 3.6]). Suppose 
/ G H„_ 24 ,.9 G H n _2j («, j < f) then 



(|M| 2j / (x) ,\\x\\ 2j g (x)) h = 5 tJ 2 2 H! (n - 2i + ^ + Ms) (/, 5 ) 



ft • 



(by the symmetry of the inner product we can assume i > j and use equation 3T ) 

■>3<n/ 



Lemma 2. Suppose f e V ( n N) and f (x) = J2i<n/2 \ \ x \\ 23 fn-ij {%) with f n _ 2 j S 



I n _2j /or eac/i j < § i/ien A^./ = 2 2j i! (n — 2i + y + /n-2j if and only if 

fn-2j = for each j > i, (that is, AJ+ 1 / = 0/ 



2 



Proof. By equation |3.l| and for each i < 

A;/(x)= £ 2 2 *(-^7-n + J + l-^-A^) |N| 2 ^ 2i / n - 2j (^). 

i<j<n/2 ^ ' * 

This is an orthogonal expansion, hence A l K f (x) equals the term in the sum with 
j = i if and only if f n -2j = for all j > i. □ 

The Lemma allows us to find the lowest-degree part of the expansion of h^ for 

/i G Vn N \ since AJ./i M = 2 l h fl ^ iei i£ i < (ii — (J® and A^/i^ = if i > /xi - /i 2 , 
by Theorem ||. We will use an elementary property of any inner product space E: 
suppose {fi : 1 < i < m} is a linearly independent set in E, for some k < to let p 
denote the orthogonal projection of E onto span{/i : 1 < i < fc}, then 

det «/<, /,))™ =1 = det ((fi, f 3 )) k i ]=1 det (((1 - p) f u (1 - p) fj))™ =k+1 . 

The proof is easy: for each i > k there are coefficients ay, for 1 < j < m so 
that p/i = i CLjifj. Let A be the to x m matrix with entries Aj-$ = Sji except 
Aji = —a,ji for I < j < k < i < m. Let G be the Gram matrix of {/i} i=1 and let 
G' = A T GA so that 

(fhfj) tor l<i<j <k 

G'a = G'a = { (fi, (1 - P) fi) = for 1 < i < k < j < to 

<(1 - p) (1 - p) £> for k + 1 < » < j < to; 
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thus det G' = (det A) 2 det G = det G and G' has a 2 x 2 block structure with in 
the off-diagonal blocks. We apply this to the projection of Pf„ onto Hf„, denoted 
by p n - We write ei for ||a;|| 2 as before. For p £ let = J2jLo^ 2 e i^Mj ( x ) 

with each h^j £ ^-2n-2i ( an d by the above discussion, h^i is a multiple of h 



for i = /ii - /i 2 ), then (I - p n ) = X^i^ 2 e i h ^,j- Tnus Ai = D n detM 
where M is the Gram matrix for |(1 — p n ) /i M : p £ Vn N \ pi > P2^- The span 
of this set is eiPj^_ 2 , an d ^2n-2 can be decomposed just like Pf„. Indeed 
|(1 - p n ) h^, : fx £ Vn N \pi = p-2 + l| is a basis for eiHf n _ 2 , since (1 - p n ) = 
ei/i M ,i, a nonzero multiple of ei/i M _ ei , for p\ = P2 + 1- Repeat the previous step 
with the projection eipn-ie^ 1 to express detM as a product. Here is a formal 
inductive argument. 

Let (A, p) be the symmetric matrix indexed by V n N ^ with entries 



M (j) (A, /x) 



(e\hx^,e\h^ t ) h for Ai - A2 = Ml - M2 = « < j 
for Ai - A 2 < min (^1 - p 2 , j) 

/Ai — \ 2 \ 

( J2 e \h\,i, e lV* ) for Ai - A 2 ,Mi - > j- 
\ l =i »=i / h 



By the symmetry M^(\,p) = if /^i — p% < min (Ai — A2, j). The matrix 
has a diagonal block decomposition (zero blocks off the diagonal) with one 

block for each set of labels |a £ vh N ^ : Ai = A 2 + i\, equivalently V^l\, for each 
i = 0, 1, . . . , n - 2, n (the set v[ N) is empty). The matrix = G, the Gram 

matrix of |/i A : A £ }• We show that det G = det by proving det = 

detM^' +1 ' for each j < n. Only the principal submatrix of AfW labeled by A 
with Ai — A 2 > j need be considered. This is the Gram matrix of a certain 
basis for e{P 2 3 I1 _ 2 ji the projection e{p„_ :( e 1 ^ maps this space onto ejH 2 3 n _ 2:) = 

span $e{h\j : Ai = A 2 + jX. Now observe that (l — e\p n -je[^\ e \h\,i = 

^ ^ * ' i=j 

A1-A2 

X) e \h\,i for A with Ai — A2 > j. The projection argument shows &etM (J > = 
i=j+i 

The determinant of the principal submatrix of M^ n ' labeled by 
|a £ vi N ^ : Ai = A 2 + i| is a multiple of D n ^. Let 

2 2l i\ \2n-2i + ^+Nn 



by Lemma Q h\.i = CiA l K h\ = 2 l Cih\^i £l for Ai = A2 + i. Further 

(e\hx,i,e{h IMti ) h = cj 1 (/i A ,i, h^) h 
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for Ai — A2 = /ii — fJ,2 = i- The correspondence a 1— > a + ie% is one-to-one from 
V { "_\ to {A e V { n N) : Ai = A 2 + »}. Thus 

D n = detG = JJ I i! I 2n- 2i+ — +iV/c) 



where d (j, iV) = #P!j N \ for j > 0. In particular, ci (1, N) = and L>i = 1. 

Theorem 3. XTie determinant D n of the Gram matrix of for the 

inner product (■,•)/, satisfies 

D n ( ( N XX d(n-i,N) 

i— 1 ^ ^ 

Proof. From the above result 

d(n-l-jM) 



— Dm 



f2n- 2i+ f +JV/«)._ 1 



In the numerator replace j by z — 1; the ratio ln N 



2n - i - 1 + f + Nk 



(2n- 2i+ f + AT/c). 
, and this completes the proof. □ 



In the ratio p n there is some cancellation by use of the decomposition Vn 
V^ N) U j/i + £1 : n e V$*\y, indeed each fi e 7^ contributes 

A*!(* + §) u m M (l") 1 



(// + £1)! (k + i) M+£i m^+e, (l w ) + 1) (k + § + Mi) # {.7 : Hj = Mi} ' 

Note that {\ N ) = N\/Y\ S>Q : Mj = s l)' so the change from (l N ) to 

m n+ei {}-) m the replacement of sl by 1! (s — 1)! where s = # {j : Mj = Mi} (except 
when /1 = (0) in which case m M (l w ) /m fJi + £l (l N ) = jf", this only affects the 
vacuous equation D± = 1). The other part of the expression for D n is a product 
with J2i=i d(n — i,N) = #v!^l\ terms (and if ju G with ^i — M2 + 1 = i then 

M G T 5 ^! + (i — 1) Ei,l < i < n). This is the simplified formula (for n > 2): 



5„= n ( a! ^( iAr )) 



-q (Mi - M2 + 1) (2n - 2 - m + na + f + Nk) 



(fix + 1) (k + I + Mi) # 0' : Mj = Mi} 



The determinant of the Gram matrix for {•,•} s is an easy consequence by Proposi- 
tion ^: 

/ / N \ \ -rf(n.iV) _ 

det((/iA,Ms)A,^<"> = ( 22 "(y + ArK ) ) ^ 
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When N > n the products in D n are over all partitions V n ,V n —i and the inner 
product formula for (h\,h^) h in Proposition || can be considered with N as an 
indeterminate because 

-C") — (,£)) 

= m„(l«">) (-1)«">(-«), W /((„)!, 

a polynomial in AT, for a € V and iV > / (a). So the formula for D n is an identity in 
the variables k, N (for N > n). It would be interesting if there were an orthogonal 
basis for so that for each element / the squared norm (f,f) h is a product of 
integral powers of factors linear in n; this is almost suggested by the nice form of 
D n but so far such a basis has not been found. The author's suspicion is that the 
problem is the lack of a sufficiently large set of commuting self-adjoint operators 
on H£. 



4. Point evaluation of harmonic polynomials 

For almost any classical orthogonal polynomial there is at least one point where 
the polynomial can be evaluated in simple terms, for example the Jacobi polynomi- 
als have such values at ±1. It is natural to ask this question about the symmetric 
polynomials h\. We will consider two special points with a very simple result for 
one and a more complicated but quite interesting analysis for the other. 

The easy evaluation is at x = (1, 0, ... , 0) because m\ (x 2 ) = if I (A) > 1; 

thus for A € V^ N) and n > 2 we have h x (x) = 2 B ( A > H) _ By T h eorem pH 

n\[K+±) N 



-l 



B (A, (n)) = (-1)" - Al n{n - A 2 - 1)! (jlL A,f 

The more interesting evaluation is at x — 1 N (one reason to single out this point 
is that the discriminant Yl%<j i x i ~ x j) vanishes there of maximum order, another 
is its role in the generalized binomial coefficients). Also for A G "Pn the value 



depends only on A (and k) but not on N, so one may as well assume N = I (A). 
Analyzing this sum as a rational function of k we observe that each term is O («~") 
as k — ► oo. We will prove that h\ (l N ) — O L(3n+i)/2j \ ^s k — ► oo 5 by expanding 
7 — — as a Laurent series at 00 and showing that the coefficients of kT 1 cancel out 

for n < i < n + [^^J (rephrasing this bound: if n = 2s then 2s < i < 3s — 1, if 
n = 2s + 1 then 2s + l<i<3s + l where s = 1,2,3,...). The cancellation is a 
consequence of the following property. 

Lemma 3. Suppose n > 2, A G Vn an d g is a symmetric polynomial in N 
variables with deg (g) < n then 

2^ — ^ — ^ (a*) = °- 
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Proof. The idea is just a slight generalization of Proposition ^. Fix j < n and let 
v G v\ be arbitrary. We consider the expansion of e™ " 'm„ in two different ways: 

e"~ J m„ = ^2 ( )"V = E B (a, fi) ( \m a 
= e™^ ^B (t,v)ui t = ^B (r, i/) m r+(n _ i ) E1 . 

bmce {m CT : a G # | is a basis it follows that X)i/cm ^ (£) = ^ ^ or an ^ ^ e 

"P^' and any v G "P.^ with < j < n. The set {/x ^ : ^ e ^j^' < j < ™f 

is a basis for the symmetric functions of degree less than n. By definition 







r w i 








r AT 




(-1) H 

I i=l 



: cr = 1/ 



Note that any terms in the sum with at least one o~j > [ii vanish, so the sum in the 
last line can be considered as a symmetric polynomial in the variables fix, . . . , hn. 
The highest degree part is clearly (n — j) \ m u (fi) . Induction shows that m T (/i) G 

span : v e vi j N) i <3< n } for each T € P,^. <i<n. □ 

The following is a rational version, with an elementary algebraic proof, of the 
asymptotic formula for T (t) /Y (t + a) for t — > oo. 

Lemma 4. There is a sequence {pj (n) : j G No} of polynomials such that pj is of 
degree 2j with leading coefficient (2 J j!) 1 /or j > and 



K n ^2p 3 { n ){- K Y 



I) 

2Jn j=0 



for \k\ > n — i, n > 0. 



Proof. The coefficients of the Laurent series certainly exist so it suffices to prove 
their polynomial property. Clearly po (n) = 1 and Pj (0) = for j > 1. There is a 
recurrence relation implicit in the equation: 

l 



1 n , 2 



OO OO s - \ 

« _n E (» + !) - Pi (»)) (-^r 3 = -" _n_1 E^ (" + !) ( " + o ) (-«)" 



2 7 

i=o j=o 7 

Extracting the coefficient of (— l) 5 k - " - - 7 in the second equation yields 

Pj (n + 1) - (n) = ( n + ^ ) Pj_i (n + 1) . 
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Arguing by induction suppose that deg (pi) — 2i for < i < j — 1 then the right 
side of the equation is of degree 2j — 1, the first difference of pj is of degree 2j — 1, 
Pj (0) = and this determines pj uniquely as a polynomial of degree 2j. Denote 
the leading coefficient of pj by Cj then the coefficient of n 2 ^ 1 in the first difference 
shows that 2jcj = Cj—x. □ 

By direct computational (n) = ^n 2 ,p2 (n) — -^n(n+ 1) (3n 2 + n — l) .p% (n) = 
n 2 (n + 1) (n + 2) (n 2 + n — l) . The explanation for the linear factors is that 
^ = (k + i + n) _ n for n, < 0, which is a polynomial of degree — n in k, hence 

Pi ( n ) = for j > — n; the recurrence is also valid for n < 0. 

Proposition 8. For A G iftere is a polynomial q\ (k) sitc/i that deg < 

Lf J — an d 



V> _ 9A («) 



n 



2; n lli=2 V." 1 r 2/Aj 

Proof. The definition of ^ implies that 

'-{H)/*= i "}"H).sH, 1 , 

thus the claimed equation is valid for some polynomial q\ such that deg (g\) < n—\i 
(since h\ \1 N ) = O (n~ n ) as k — > oo). Assume \k\ > n (only a temporary restriction 
since we are dealing with an identity of rational functions) and expand each term 
in the sum for h\ (l N ) by Lemma [| to obtain 

_. I W 

For any s > the sum ^ (— k) J} p Qi (//j) is a symmetric polynomial 

a£N™,\a\=s i=l 

of degree 2s in /i. In the Laurent series /i A (l ) = J2 < jLn a 3 K ~ : ' ^ ne coefficients 



a 



J r- 

for 2 (J — n) < n by Lemma that is a,j — for j < |n and ft. A (l^J 
O ( re -L(3n+i)/2jy Thusdeg(g A )-(2n- Ax) < - (n + L^J) and deg (g A ) < [f J - 
A,. " □ 

The next step is to find an explicit form for q\ by imposing [§J — Ai+1 conditions 
(the number of coefficients); these are obtained by a calculation exploiting the fact 
that h\ (l N ) has simple poles at k = h—j for Ai + 1 < j < n. There are more points 
than needed so we consider only the residues at n = | — j for [^y^-J + Ai < j < n 
(note that [^J + L 1 ^] = n so there are [§| — Ai + 1 points in this list). Only the 

P art E \ j7 ^ T\ : A ^ M, Mi > L^J + A i f of t hc sum for ^a (l^) has poles 
at these points. 

Suppose 1 < j < s then the residue of -; p— at K = \—j equals ■ 



U-iy.(s-jy. 

(remove the factor n — h + j and substitute k = i — j). 
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Theorem 4. For A e vi N) let k (A) = Ai + [^J then 



hx (l N ) 



2 

2 -n 



v b(x^) * niI 2 (i-j+^) Ai _^ (-ir fc < A ) 

A.J^A) ^ ! A) (^- fc (A))! 
Proof. By Proposition ^ 

^ = 7^ — ^ ^fr? 

l K + 2J A+S£1 j=fe(A) K 

with s = k (A) — Ai — 1 for certain coefficients {r^}. Multiply both sides of the 
equation by (k + |) A+S£ to obtain 

V 5 _ 2 -»v i3(A ^ ) ^ + D fc( A)- 1 A (^ + DA, 

Note that = ( K + | + Mi) a ( for * > 2 and A ^ /i). For j > k (A) 

the term corresponding to some /i has a pole at k = \ — j exactly when j < /ii 

-B (A f— lV — * 

(and A ^ /z); the residue of this term is } — - — — (f — x 

M- U - !) ! (Mi - J)' 1 ; 

AT 

n (i-i+w)^-^- □ 

4=2 

Technically this is not a summation but a transformation formula, but it does 
use fewer terms than the original form of h\ (l N )- If Ai is close to S then the 
formula uses relatively few terms. There is only one term for A = (1,1) or (1,1,1) 
for n = 21 or 21 + 1 respectively. It is easy to show B ((I, I) , (21)) = 2 (—if and 
B ((1, 1,1), (21 + 1)) = (-1)' -1 (2Z + 1) . Then the theorem implies h {u) (l N ) = 

( 22 ^ ! (- + I)(2M)) _1 and Wi") = ( 22i C-^K^+Dc^w))" 1 - 

value (l^) can also be found directly by using Kummer's terminating sum for 
2 F 1 (-21, b;l-2l-b; -I) with b=\-K-2l. 

5. Concluding remarks 

There is an isometry between the space of -invariant polynomials with the 
inner product {■,-)$ and polynomials with the L 2 -inner product for the measure 

(r (f + Nk) /r (k + i) W ) nilj 2/r 1/2 d?/i ■ ■ ■ <*ttv-i, on the simplex 
jy G : J2iLi Vi = ^-iVi — eacn induced by the correspondence y = 
(xf,... t^jv)- Then is isomorphic to the space of polynomials of degree n 
orthogonal to all polynomials of lower degree. The basis ih\ : A € "PrV^ j maps to 

a basis for the polynomials symmetric in (y\, ... ,tjn) ■ 

If one gives up the -Bjv-invariance then there is no problem in constructing a nice 
orthogonal basis for H n . This basis consists of products of Jacobi polynomials; a 
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conceptual derivation in terms of simultaneous eigenfunctions of a set of commuting 
self-adjoint operators can be found in Theorem 2.8]. 

In this paper we constructed a basis for B^v-invariant spherical harmonics by 
introducing a new basis for the symmetric functions. The determinant of the Gram 
matrix of the basis was explicitly evaluated. Long ago it was found that the wave 
functions of electrons in a crystal with cubical symmetry have (energy level) degen- 
eracies for n > 12; in our notation vf ] = {(1, 1)} , {(1, 1, 1)} , {(2, 2)} , {(2, 2, 1)} for 
i = 2, 3, 4, 5 respectively but Pg 3 "* = {(2, 2, 2) , (3, 3)}. We have not found a natural 
and constructive way of orthogonally decomposing when dimH^ > 1; a candi- 
date for such a decomposition is the self-adjoint operator J2i<i<j<N { x i^j ~ XjT>i) 
but its eigenvalues are irrational over the field Q (k) of rational functions in k. Yet 
the methods in this paper should give some insight into the problem of constructing 
invariant harmonics for the general -B-type spin Calogero-Moser model. 
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